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A formulation of Einstein equations is presented that could yield advantages in the study of 
collisions of binary compact objects during regimes between linear-nonlinear transitions. The key 
idea behind this formulation is a separation of the dynamical variables into i) a fixed conformal 
3-geometry, ii) a conformal factor possessing nonlinear dynamics and iii) transverse-traceless per- 
turbations of the conformal 3-geometry. 
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In the two-body problem in General Relativity, involving black holes and/or neutron stars, one can identify three 
main phases according to the extent to which nonlinear effects play a fundamental role. When the two bodies are 
sufficiently far apart, the binary system is in an Inspiral phase. During this phase, a Post-Newtonian approach to 
gravity provides a suitable framework. Loses due to gravitational radiation drive the binary system to the next stage, 
)—} , the Merging phase. During the merging, which includes the last few orbits and coalescence of the binary system, 
■ nonlinear effects become essential and eventually dominate the dynamics of the system; thus, the full set of Einstein 
fields equations is required. The last stage of the problem is the Ringdown phase. Here once again, the two-body 
y—i . problem can be treated perturbatively. In this situation, however, the system can be investigated as perturbations 
about the outcome of the collision, namely a single black hole or neutron star. This point of view is also known as 
. the "close-limit" and has been shown to be a remarkable good approximation ([J . 

From the computational point of view, it is safe to say that the machinery, necessary to study the Inspiral and 
Ringdown phases of the problem, is available and reasonably well understood. On the other hand, in spite of advances 
in computer technology and numerical algorithms, the Merging phase of black hole and neutron star collisions remains 
^\ I essentially unsolved; that is, simulations that bridge the "gap" between the Inspiral and Ringdown stages are not yet 
ON ■ available. The expectation is that this situation will prevail for the next few years. Furthermore, even if full nonlinear 
' numerical simulations of the collision of compact objects are achieved they will most likely not be able to cover, at 
^j*. least initially, the dynamical time-scales necessary to grab initial data from the Inspiral phase, evolve it through the 
Merging phase and provide initial data for Ringdown calculations. To get a sense of the formidable task involved, one 
only needs to keep in mind that, for numerical approaches using asymptotically inertial coordinates, the dynamical 
I* ■ time-scale is ~ M. This time-scale should be compared with ~ 1200 M, the time-scale of numerical evolutions to 
. £^ | cover the Inspiral- Ring down gap. Because of this foreseen limitation, it is important to develop methodologies that 
^ . directly address the intermediate regimes between the main phases discussed above. 

The Inspiral- Merging intermediate regime has recently received significant attention. In the context of black hole 
collisions, this is called the "intermediate binary black holes" problem by Thorne and collaborators ||. They have 
proposed an approach in which numerical calculations are performed in co-rotating coordinates with the binary system. 
In this co-rotating frame, the metric evolves on long time-scales, namely those of the inspiral scales. Furthermore, 
gauge conditions are chosen such that gravitational degrees of freedom that are not excited by radiation remain frozen. 

Finding a reduced formulation of Einstein equations for the Inspiral- Merging as well as Merging- Ringdown inter- 
mediate regimes is the focus of this paper, namely a formulation tailored to the transitions between full non-linear 
and completely linear phases. Our working premise will be that, during those linear-nonlinear (LNL) transitions, 
"democracy" does not apply when determining the framework needed to handle the gravitational fields. That is, we 
expect that nonlinear effects are not equally turned on and off for all the field variables. There is evidence for such 
behavior. In the numerical study of neutron star binaries, Wilson and collaborators Q|J used a simplified version of 
the Einstein equations in which only a reduced set of metric variables have nonlinear dynamics, with the remaining 
fields fixed in time. Although, a debate arose on one of the predictions from this work, namely the collapse of the 
neutron stars before their coalescence (see for a resolution to this debate), this approach provides a reasonable 
approximation in the regime when the dynamics is dominated by volume (conformal) deformations and the energy 
radiated per orbit is much smaller than the energy of orbital motion. Our LNL formulation of Einstein equations 
builds upon Wilson's idea, allowing for more general background space-times and the presence of linear perturbations. 
The LNL is constructed by identifying, among the dynamical variables, those quantities that can be kept fixed, those 
that can be treated perturbatively and finally those that follow nonlinear dynamics. The hope is to arrive to a reduced 
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system of Einstein equations that is not only more amenable to be handled numerically but also allows one to gain a 
better insight of the problem. 

As customary in Numerical Relativity, the starting point is the 3+1 or ADM form of Einstein equations |F3]. Under 
the 3+1 scheme, the space-time is viewed as the time history of a foliation of space-like hypersurfaces with intrinsic 
metric h ab , with Latin and Greek letters denoting spatial and space-time indices, respectively. Since most of our 
expressions will eventually involve conformal quantities, we will reverse the usual convention and denote physical 
quantities with hats and conformal quantities without hats. The 3+1 line element takes the form: 

ds 2 = -a 2 dt 2 + h ab (dx a + (3 a dt)(dx b + (3 h dt) , (1) 

with d the lapse function characterizing proper time along the normals to the hypersurface, and a the shift vector, 
representing the freedom of relabeling spatial coordinate points in that surface. With the space-time metric given by 
(Jl|), Einstein equations take the form: 

R + K 2 - k ab k ab = (2) 
V b K ab - V a K = (3) 
B h ab = -2d K ab (4) 
d K ab + 2 d K ac K c b - a K ab K — a (R ab - d ab - d a d b ) (5) 

Above, R ab is the Ricci tensor constructed from the 3-metric h ab , R = R a ° is the scalar curvature, and K ab denotes 
the extrinsic curvature, with K its trace. Also, a a — V a In d is the acceleration of Eulerian observers, and we introduce 
the tensor a ab = V a d b — a ba , with a = d a a = V a a a . Covariant differentiation with respect to h ab is denoted by 
V , and d = d t — £g. The operator d -1 o is the derivative with respect to proper time along the normal to the 

space-like hypersurfaces, and V L = d + (3^ is the time vector tangent to the congruence of world lines of coordinate 
observers. Equations (0) and (^) are respectively the Hamiltonian and momentum constraints. Equation (Q), on the 
other hand, follows directly from the definition of the extrinsic curvature K ab . For simplicity, we are only considering 
vacuum space-times; however, it is straightforward to include matter sources. 

The first step to derive the LNL formulation of Einstein equations involves splitting the evolution equations for the 
metric h ab and extrinsic curvature k ab into trace and trace-free equations. That is, using K ab = A ab + h ab K /3 with 
A a a — 0, Eqs. (^]|) are rewritten as: 

R + \k 2 - A ab A ab = o (6) 

V b A ab - \\/ a K = (7) 

d \nh 1/2 = -aK (8) 
B k-ak 2 = a(R- a- d a d a ) (9) 
h 1 ^d (h- 1 / 3 h ab ) = -2aA ab (10) 

d A ab + 2 d A ac A c (,--q A ab K = a (R ab - d ab - a a d b ) (11) 

where angle brackets denote traceless symmetrization, i.e. (T ab ) = T^ ab ^ — | h ab T with T = T a a . 

In numerical relativity, conformal transformations have been extremely useful, in particular in the construction of 
initial data || to single out those "pieces" within h ab and A ab that are fixed by the constraints (JsJ) and (]]). Recently, 
conformal transformations have been considered also in formulations of the evolution equations &P~11,- There are 
some indications that these formulations exhibit improvements in numerical stability properties p-|ll|. Here, we 
introduce conformal transformation as a tool for factoring out nonlinearitics. There is, however, a certain degree of 
ambiguity regarding the particular form of these transformation rules. Once the conformal transformation rule for the 
spatial metric h ab is fixed, there is no specific recipe for choosing the remaining transformation rules for the extrinsic 
curvature K ab , lapse function d and shift vector (3 a . For instance, in the construction of initial data, York's choices 
H were mostly aimed at simplifying the task of solving the constraints (|J) and (Q). After several attempts, we found 
that the most convenient conformal transformations that yield a LNL factorization are: 

h ab = e 4n h ab (12) 
/?/ 2 = e 6n (13) 
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K = e~ 2n K 



a 

/3 a = (3 a 



2H 

e a 



(14) 
(15) 
(16) 
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Conformal transformations similar to the above were used by Geroch jlj] to study the structure of the gravitational 
field at infinity. The conformal transformation ([l2]) for the spatial metric is just the standard rule used in most of the 
literature. The transformation (|l^) for the determinant of h a b implies that h = 1 . This choice is always possible 

and simplifies the form of the e qua tions while keeping the procedure general. A consequence of h = 1 , together with 
the conformal transformations ( |12| - |l7| ) , is that h a b is a tensor density of weight —2/3 and A ab a tensor density of 
weight —1/3. Similarly, if is a scalar density of weight 1/3 and a a scalar density of weight —1/3. Notice also that, 
contrary to the common practice of not assigning a conformal transformation to the trace of the extrinsic curvature 
K, we have chosen a conformal transformation for K such that K ab — e 2n K ab . It is also important to point out that 
with the above conformal transformations, the line element (Q) takes the form 



ds z 



„4 fir 



-a 2 dt 2 + h ab (dx a + f3 a dt)(dx b + f3 b dt)] 



(18) 



Therefore, the space-time metric has the same transformation rule as the spatial metric. Furthermore, the transfor- 
mations ([L2-17) yield quantities a and j3 that are the lapse function and shift vectors of the conformal space-time, 
which would not necessarily be the case in general. 

Given (Ojl^), the following conformal relationships are also obtained: 



-Rah = Rab - 2 UJab + 4 UJ a UJ b - 2 ll ab (uJ + 2 UJ C LU C ) 

R = e- 4n (R-8uj-8 uj a uj a ) 
a a = a a + 2 uj a 

dab = a a b + 2 Lu ab - 4 0( a u b ) -8w a w 6 + 2 h ab (a c uj c + 2 u c lu c ) 
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(a + 2 u + 2 a a uj a + 4 uj a u a ) 



(19) 
(20) 
(21) 
(22) 
(23) 



where w a — V fi and oj ab = V a u b — ui ba with to = u> a a = V a Lo a . We also have that d — d and (V°/3 fc ) = 
e -4fi (y a f3 b y With the set of conformal transformations (|l2-23), Einstein equations (pHTl]) can now be rewritten as: 



R+-K 2 - A ab A ab 



8(w + Lu a u a ) 
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where 



V b A ab - - V a K = -Auj b (A ab + - h ab K) 



d o n = --aK 
6 

2 

d K --aK 2 = a(R-a- a a a a - n) 
3 

d h a b = ~2aA ab 

2 

d Q A ab + 2a A ac A c b - - a A ab K = a {R ab - a ab - a a a b - U ab ) 



U ab = 4 uj ab - 8 uj a uj b + 2 h a b (w + 4[i) c w c + a c uj c ) 
II = ^ afc n Qb = 10w + \Quj a u a + 6a a u a . 



(24) 

(25) 

(26) 

(27) 
(28) 
(29) 



(30) 
(31) 



In the form given by Eqs. (26-29J), the Einstein evolutions equations are split into two sectors: (CI, K) and (h abl A ab ). 
The (fi, K) sector is kept fully nonlinear, and, as we shall see later, approximations are introduced only when hand ling 
(h abl A ab ), specifically in the transverse components of those fields. We start by rewriting the evolution equation ( |2S| ) 
of the conformal metric as 



d t h ab - £ph ab = d t hah - 2 (V„ j3 b ) 



-2aA ab 



(32) 



where the substitution £ph ab = 2 (V a Pb) in the second equality uses the condition that h a b is a tensor density of 
weight —2/3. Next, we decompose the tensor aA ab into its transverse and longitudinal parts, namely 
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a A ab = T ab + L ab , (33) 

such that V b T ab = T = and L ab = (V a £ fc ). Substitution of @ into (||) yields 

d t h ab = 2 (V„ (A - 6)} - 2 T ab . (34) 

Since our goal is to treat the transverse-traceless fields perturbatively, it is then natural to choose the shift vector (3 a 
such that it eliminates the longitudinal terms in (|32"|); that is, (3 a — £ a . With this choice, Eq. ( |3~4| ) becomes 

d t h ab = -2T ab (35) 
iab = (V a /3 6 ). (36) 

Notice that this choice of shift vector implies that the conformal space satisfies a gauge condition belonging to the 
class of "radiation gauges" , namely 

V b d t (h N h ab ) = 0, (37) 



with h = 1 and N = 1/3. Condition (37), however, does not translate to a radiation gauge condition in the physical 
space-time since 

V b d t (h 1 / 3 h ab ) = V b d t h ab + Wuj b d t h ab , (38) 

and in general the second term of this equation does not vanish. 

From the decomposition ( |33|) and the transverseness of T ab , it follows that 

V b L ab = \7 b (aA ab ), (39) 

which can be rewritten using Eq. (Eq) as 



y b L ab - (a b - 4 uj b ) (L ab + T ab ) = - a {V a K -2u a K). (40) 

o 

This equation is simply the momentum constraint and must be satisfied independent of the choice of shift vector. 
However, when the specific gauge choice ( p^ ) is made, Eq. ( (40| ) is alse viewed as a gau ge condition from which the 
shift vector is constructed. A more geometrical interpretation of the gauge condition ([l0|) is obtained by noticing that 
this condition, also known as "minimal stria n g auge," can be derived from a global minimization of the time rate of 
change of the metric in the conformal space |13j, that is, from the variational principle 

SS[a, l3 a ] = S [ d t h ab d t h cd h ac h bd dv = 0, (41) 



varying (3 a . We must emphasize that condition (39), or equivalently (|40|), is a minimal strain gauge not applied to 
the physical space-time but to the conformal foliation. 

The freedom of imposing a gauge condition in connection with the lapse function a still remains. For instance, 



if the choice is to have maximal conformal slices, i.e. K = 0, Eq. (27) becomes an elliptic equation for a, namely 
a + a a a a = R — II. This condition also implies from Eq. ( |26| ) that d o 0, — 0; that is, the conformal factor £1 is 
constant along the normals to the hypersurfaces. Other possible choices are the use the generalized harmonic gauge, 
a a + a K — a /, with / an arbitrary function, or a lapse function constructed from the variational principle ( f4l| ) by 
varying a instead of (3 a . The latter approach yields 

aA ab A ab = (V a (3 b )A ab , (42) 



which with the help of the Hamiltonian constraint (24) can be rewritten as: 



L ab (T« b + L« b ) 

a = t — o ; r ■ (43) 

R + | K 2 - 8 (lo + uj a oj a ) K ' 

Notice that with this choice of gauge, the equations for the lapse and shift become a coupled set of equations ||. 

So far, no approximations have been made. The theory remains completely nonlinear. The only condition is that the 
shift vector satisfies Eq. ( f40| ) , a condition introduced to eliminate the longitudinal terms from the evolution equation of 
the conformal metric. We now assume that, in the physical systems where our LNL form of Einstein equations apply, 
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the non-linearities are dominated by conformal (volume) deformations. As mentioned before, numerical simulations of 
neutron star binaries [Q,[5 14 , lj| have shown that there is indeed a phase during the inspiral where this approximation 
holds. The main difference in the present work from those studies is that we allow for the presence of perturbations 
entering via shear deformations of a conformal background. That is, we decompose the conformal metric as 

K b = ^ab + *afc , (44) 

with dt^ab = and $ a b a perturbation. An immediate consequence of the assumption (^J) is that from Eq. ( |j5|) 

d t y ab = -2T ab , (45) 

so T a b becomes a perturbative quantity, in agreement with our goal of treating the "radiative," transverse-traceless 
fields as perturbations. Furthermore, having dtjab — and T a b transverse-traceless implies that Eq. (Eft) yields 



d t ^ = 

av ft * Qb = o. 



(46) 
(47) 



where "J = 7 afcv E , a 6- Above and from now on, tensors have their indices lowered and raised with the background metric 
7ab, and tildes denote quantities in connection with this metric. A direct implication of Eqs. (^6|) and ([47j ) is that, 
if we choose initially the perturbation of the metric to be transverse and traceless, these conditions are preserved by 
the evolution equation (^). We will take advantage of this property and assume that initially Vb^E ,ab = '5 = since 
this choice simplifies the equations significantly. 

At this point, we have the Hamiltonian constraint (f24|), the gauge condition (|io| ) on the shift vector and the nonlinear 
evolution equations (|2^) and (|27|) for the conformal factor and the trace of the extrinsic curvature. In addition, we 
have the evolution equation (f45|) for the perturbation to the conformal metric; therefore, what remains is to obtain 
an evolution equation for T ab . 



From Eq. (|45j), we have that 



where we have defined 



Oo^ab 



-2 T, 



P^ab 



(48) 



Tab — T a b + - £p^ab ■ 



On the other hand, Eq. (33) can be rewritten as: 

a A a b = L a b + T a b 



V Q (3 b ) + T a b 
£ph a b + Tab — £plab " 



L a b 
1 

(V a Pb) + T a b = Lab + %b ■ 



(49) 



1 



£^ab + Tab 



(50) 



where we have used that "/ a b is a tensor density of weight —2/3. Substitution of fl50| ) into (29) then yields 

d Q (a A ab ) - (a Q - 4 d o n) (a A ab ) + 2 (a A ac ) (aA c b ) = 
d L a b - (a„ — 4 d fl) L a b + 2 L ac L c b + 

d Tab - (a a - 4 d o n) Tab + 2 {L ac T c b + L bc T c a ) - 2 L ac L db * cd = 

Oi 2 {R a b ~ dab — a a df, — ±I t,} , 

where 



(51) 



(a a b) 
(a a ab) 
(R a b) 

(Rab) 



(aab) + - {lab^ cd - l cd ^>ab)a cd - S c a b a c 
(a a a b ) + ^(7afc* c<i - l cd ^ab)a c h d 



(Rab) + ^(lab^ 



l Cd ^ab)Rcd ~ g V c V C * ab + R c{a ^ b ) ' 



Rr.ahd^ C 



(U ab ) + g (7afc* - 7 *ab)n c(i - 4 S c ab Cj c + 2 V ab (u + 4 Cj c Cj c + h c u c ) . 



(52) 
(53) 
(54) 
(55) 
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with 



S a bc = g (V 6 * a c + V c * a b - V a * 6c ) ; 
therefore, to zero-order Eq. ( |5l|) reads 

9 L ah - (a G - 4 9 f2) L ah + 2 L QC i c 6 = a 2 (R ab - a ab - a a a b - U ab ) 



(56) 



(57) 



These evolution equations are redundant because, as mentioned before, the gauge conditions (|4^) that determine L ab 
are also the momentum contraints, which in turn preserve Eqs. j57j). To first-order in the perturbations Eq. j5l| ) 
yields the evolution equation for the perturbation T ab : 

d T ab -(do -4 o n) T ab + 2 {L ac T C b + L bc T C a )~2 Laddb ^ = 

- ^ VcV^ah + R c[a ^b) c + Rcabd^ cd + S c afc (5 c + 4 Cj c ) - 2 * ab (£ + 4 U) c Cj c + h c Cj c ) 



+ g(7ab* - 7 ^ab){Rcd ~ a cd -h c h d - n cd ) 



By combining this equation with Eq. (|48|), one obtains 



where 



— d -d + V c V c )V ab = J ab , 
a a 



^a 2 J ab = ^aKTab~2 (L ac T c b + L bc T c a ) + 2 L ac L db * 



(58) 



(59) 



cd 



Rc(a^b) c + Rcabd^ cd + S c ab {a c + 4w c ) - 2 <f afc (w + A£o c Cj c + h c Q c ) 

+ |(7a6* Cd - l Cd ^ab){Rcd ~ a cd ~h c h d - ft c( 



(60) 



That is, the perturbation $ a b obeys a wave equation in the conformal space-time. 

One more simplification of the LNL form of Einstein's equations is possible: neglecting the back-reaction of the 
perturbative fields (^> a b,7~ a b) on the nonlinear (fl,K) and gauge (a,/3 a ) fields. Under this approximation, covariant 
differentiation and the raising and lowering of indices are performed with the background metric "f a b- 

In summary, the following picture emerges. Starting from the standard 3+1 formulation of Einstein equations in 
terms of variables (h a b, K ab ), we have constructed a formulation in which the field variables are divided into: 



1. Nonlinear fields (Q, K) with their dynamics given by 

d fl = --aK 
6 



d K--aK 2 =a(R~a- a a a a - ft) 
3 



(61) 
(62) 



2. Perturbative fields (^ a bi Tab) with equations of motion 

d ^ ab = -2Tab 
1 



d T ab - a T ab = -^a 2 (V c V c # & - Jab) 



(63) 
(64) 



3. A background conformal metric j a b such that dtj a b = 0. 
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4. And a shift vector f3 a satisfying the condition 

V b L ab - (o 6 - 4 u b ) L ab = \ a (V a K -2Q a K). (65) 

The conditions and assumptions used to derive this LNL system of equations are: 

1. The field variables are conformally transformed according to ( jl2[]l7| ). 

2. The shift vector is chosen to provide the longitudinal part of aA ao , which in turns implies that it must satisfy 
Eq. @. 

3. The conformal metric h a b is decomposed into a background metric 706 plus a transverse-traceless perturbation 

4. The back-reaction from the perturbations is ignored. 

As mentioned before, the only gauge freedom remaining is that in connection with the lapse function a. Furthermore, 
if a fully constrained evolution is preferred, the conformal factor could be obtained from the Hamiltonian constraint 
(p4|). When in addition one assumes that the metric perturbation 'fab, and consequently Tab, could be neglected, the 
LNL system reduces to a nonlinear system similar to that used by Wilson and collaborators [Q^) for studying neutron 
star collisions. The differences would be due to our particular choice of conformal transformations. 

In conclusion, with the help of suitable conformal transformations and a shift vector gauge condition, we have 
constructed an approximate formulation of Einstein equations that could yield advantages in the study of gravitational 
systems where radiative, transverse-traceless variables can be treated as perturbations of a fixed, background spatial 
metric, with the remaining fields possessing nonlinear dynamics. Since this separation depends on the choice of 
conformal transformations, it remains to be tested whether the conformal transformations used in our derivation of 
the LNL system of equations indeed constitute an appropriate approximation. When the LNL system is applied to the 
Inspiral- Merging intermediate regime, an interesting approach suggested in Ref. [gj that could place our approximation 
in a more firmer grounds, is to perform the calculations in a co-rotating frame. The objective under this approach is to 
find a spatial vector 6 M , i.e. 6 M n M = with n M the normal to the constant t hypersurfaces, such that £ M = a + (3^ + 
is an "almost Killing vector field," namely £^g flv ~ 0. 
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